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ABSTRACT 



We present the generalization of a recently introduced modified gravitational potential for self-gravitating fluids. The use of this poten- 
tial allows for an accurate approximation of general relativistic effects in an otherwise Newtonian hydrodynamics code also in cases 
of rapid rotation. We test this approach in numerical simulations of astrophysical scenarios related to compact stars, like supernova 
core collapse with both a simplified and detailed microphysical description of matter, and rotating neutron stars in equilibrium. We 
assess the quality of the new potential, and demonstrate that it provides a significant improvement compared to previous formulations 
for such potentials. Newtonian simulations of compact objects employing such an effective relativistic potential predict inaccurate 
pulsation frequencies despite the excellent agreement of the collapse dynamics and structure of the compact objects with general 
relativistic results. We analyze and discuss the reason for this behavior. 
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1. Introduction 

Rapid rotation plays an important role in numerical simula- 
tions of a number of astrophysical events connected to stel- 
lar core collapse and neutron star evolution, such as accretion- 
induced collapse (IDessart et al.ll2006ll2007al). the coUapsar sce- 
nario for long gamm a-ray bursts (MacFady en & Woo slev 1999; 
IDessart et al.l l2007bh . the phase-transiti on-induced coll apse of 
a rotating neutron star to a quark star dLin et al.ll2006l) . or ro- 
tational instabilities like the low instability in a col- 
lapsin g stellar core (ICentrella et all l200lh. proto-neutron star 



(PNS: IShibata & Sekiguchil 120051: lOtt etall I2005L l2007al). o r 
cold neutron star dShibata et alJ 120021. 120031: lOu et alJ 12004 



lOu & Tohlindl200"6l NS;). In all of these cases, as compact ob- 
jects are involved, general relativistic (GR) effects should be 
taken into account properly. Moreover, extremely high radial 
and/or rotational velocities may be encountered, and thus a rela- 
tivistic description of self-gravity and of the flow dynamics may 
be required. However, as an elaborate treatment of the micro- 
physics in the form of neutrino transport and a microphysical 
equation of state (EoS) is crucial for some of the scenarios, and 
as multi-dimensional GR codes for numerical simulations may 
either be unavailable (as in the case of detailed Boltzmann neu- 
trino transport) or computationally very expensive, most previ- 
ous investigations relied on the Newtonian approximation for 
the strong gravitational field of the compact objects, and on 
Newtonian flow dynamics. 

A straightforward and easy to implement way of overcoming 
the first, most serious deficiency has been suggested recently for 
supern ova core collapse involving at most moder ately rapid ro- 
tation dRampp & Jankal 120021: iMarek et al.ll2006l) . The strategy 
presented in these studies is to treat relativistic gravity by us- 
ing an effective relativistic gravitational potential, while retain- 
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ing the equations of Newtonian hydrodynamics. This allows for 
a simple upgrade of existing Newtonian hydrodynamics codes 
to treat relativistic gravity in an approximate but quite accurate 
manner. Extensive tests of this approach have been presented 
for simulations of supernova core collapse in spherical symme- 
try in cluding Boltzmann neutrino transport and a microphysical 
EoS dLiebendorfer et all 120051 iMarek et al.ll2006l) . The results 
of these tests show an excellent agreement with those obtained 
by corresponding fully relativistic ID radia tion hydrodynamic 
simulations. In addition, Ma rek etall (2006) (Paper I hereafter) 
who studied rotational core collapse with a simplified EoS find 
good agreement between the results obtained with an axisym- 
metric Newtonian hydrodynamics code extended by an effective 
relativistic gravitational potential and those obtained with a GR 
code for cores with slow and moderate rotation. However, for 
rapidly rotating cores the eff ective relativistic potentials consid- 
ered by Mare keTaT] d2006l) cannot reproduce the dynamics of 
the collapse correctly, and e.g. underestimate the maximum den- 
sity at bounce by more than 30%. In particular, their potential A 
yields excellent results in the spherical limit, but clearly fails for 
rapid rotation, while potential R performs slightly better in the 
latter case, but persistently overestimates the maximum density 
during the evolution for slow or no rotation. 

The results of Paper I suggest that the concept of an effec- 
tive relativistic gravitational potential can in principle be used 
to model compact stars with rapid rotation. However, additional 
corrections to account for stronger centrifugal forces are neces- 
sary. To this end, we introduce here a new effective relativistic 
potential which contains correction terms for the case of rapid 
rotation. We assess the applicability of the new potential to the 
scenarios listed above, by presenting comparisons with both GR 
simulations and Newtonian simulations using two of the effec- 
tive relativistic potentials described in Paper I for a range of core 
collapse and NS test models. We explore the new effective rela- 
tivistic potential in the context of stellar core collapse with non- 
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negligible rotation using either a simplified EoS or a microphys- 
ical equation in combination with an approximate deleptoniza- 
tion scheme. Furthermore, in order to test whether the new effec- 
tive relativistic potential is also appropriate for compact equilib- 
rium configurations, we consider rotating equilibrium polytropes 
which serve as models for a cold NS. We investigate pulsations 
of these models and briefly discuss the applicability of effective 
relativistic potentials to NS asteroseismology. 

All of the above mentioned astrophysical events are inter- 
esting sources of gravitational radiation. Therefore, the extrac- 
tion of the gravitational wave (GW) signal from numerical sim- 
ulations (possibly using effective relativistic potentials) and the 
prediction of wave templates for current and future GW de- 
tectors of interferometric and resonant typ e is an important is- 
sue. T o this end, we discuss the results of Obergaulinge r et all 
(2006), who used one of the effective relativistic potentials from 
Paper I and report an overestimation of the GW amplitudes of 
up to 50% compared to simulations in GR even for slow rota- 
tors. We resolve this discrepancy by introducing a formulation 
of the Newtonian quadrupole formula that is better suited for 
both Newtonian simulations with an effective relativistic poten- 
tial and fully relativistic simulations. 

Throughout the article we use geometrized units with c = 
G = 1. 



2. Effective relativistic potential 



and the spherical Newtonian mass as 



In previous work ([Rampp & Janka 2002, and Paper I) it was 
proposed to approximate the effects of GR in a Newtonian hy- 
drodynamics code by solely substituting the Newtonian poten- 
tial by an effective relativistic gravitational potential <t> e ff which 
mimics the deeper potential well of relativistic gravity. We 
first recapitulate some properties of the Newtonian potential in 
Se ct. 12. 11 of the effective relativistic (TOV) potential introduced 
by [Rampp & Janka d2002l) in Sect. 12.21 and of its most accu- 
rate improvement for slow and moderate rotation from Paper I in 
Sect. 12. 31 Then we present in Sect. |2.4| the new effective relativis- 
tic potential including rotational corrections for rapidly rotating 
configurations. Finally, in Sect. |2.5| we summarize the treatment 
of GR in the conformally flat approximation, which is applied 
both in some previous studies and also in the present investi- 
gation of effective relativistic potentials for supernova core col- 
lapse for the purpose of comparison. 

2.1. Newtonian gravitational potential for a self-gravitating 
fluid 

In a genuine Newtonian simulation, which we refer to as Case N 
in the following (see also Paper I), the Newtonian hydrodynamic 
equations are solved in conjunction with the (multi-dimensional) 
Newtonian gravitational potential 



<D N = - f dr' d0' dip' r' 2 sin & — - — 
Jv \r-r'\ 



(1) 



where p is the (rest-mass) density of the self-gravitating fluid. 
This potential is the solution of the Poisson equation 



AO N = Anp. 



(2) 



For later use we also define the spherical part of the 
Newtonian potential by 



f°° dr' _ 
%rM = - — m N 

Jr f 



f dr' , 

Jo 



m^(r) — An I dr r (p) , 



(4) 



where any hydrodynamic quantity enclosed by brackets, (...}, 
is the angular average of its multi-dimensional counterpart. 

2.2. TOV potential for a self-gravitating fluid 

For devising an effective relativistic potential for a self- 
gravitating perfect fluid we require that for a spherically 
symmetric configuration it should reproduce the Tolman- 
Oppenheimer-Volkoff (TOV) solution of hydrostatic equilib- 
rium in GR. Then a first approximation to include relativistic 
effects in a Newtonian code is to replace the spherical part O(r) 
of the gravitational potential 0(r, 6, (f>) by the spherical TOV po- 
tential T ov (Ram pp & Jankall2002h : 



O r = O n - <E> N + TO v- 



(5) 



As in Paper I we refer to this effective relativistic potential as 
Case R (or "reference" case). The spherical TOV potential is 
defined as 



X°° dr' (h) 
— (Hjov + 4tit' 3 (P)j -j—, 
r-rov 



(6) 



where h = 1 + e + P/p is the relativistic specific enthalpy, e is the 
specific internal energy, and P is the fluid pressure. The spherical 
TOV mass is given by 



mxov(r) = An f dr' r' 2 (p) (1 + <e» , 
Jo 

and 



TtovM = \ 1 + <Vr> 



2mxov 



(7) 



(8) 



with v,- being the Eulerian radial fluid velocity. 



2.3. Improved potential for a self-gravitating fluid without 
rotational corrections 



Mare k et aD (120061) found in their study that the effective 
relativistic potential O r overestimates the GR effects due 
to the straightforward combination of relativistic gravity and 
Newtonian dynamics. Thus, they proposed another effective rel- 
ativistic potential 



O a = O n - O n + Otova, 



(9) 



labeled as Case A, which involves the modified spherical TOV 
potential 

T ov,A(r) = -J ^ T (m TOV ,A+47rr' 3 <P))^-. (10) 



TOV.A 



The latter differs from OtovM in Eq. (O by the modified spher- 
ical TOV mass 



(3) m T ov,AM = 4tt 



f dr' r' 2 
Jo 



<p)(l +<e»r TOV ,A> 



(11) 
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which compared to nijoy(r) in Eq. (0 contains an additional 
factor 



TOV.A 



(r) 



1 + <v r > 2 



2m TO v,A 



(12) 



Note that this factor also appears explicitly in the definition of 
the modified spherical TOV potential Otov.aM m Eq. ( fTOb - 

Since F T ov,a < 1 the modified spherical TOV mass 7m T ov,a is 
smaller than m-rov- The fact that the right hand side of Eq. (fTTT i. 
which defines the modified spherical TOV mass m-rov,A, depends 
itself on mrov.A' because of the factor Fjov.a, causes no problem, 
as mxov.A can be computed in a straightforward way by means 
of an ODE integration. 

As demonstrated in Paper I, for simulations of rotational stel- 
lar core collapse the usage of the effective relativistic potential 
® a results in collapse dynamics and a PNS structure that are 
very close to the outcome in GR, at least in the regime of slow 
and moderately fast rotation. 

2.4. Effective relativistic potential with rotational corrections 

Following the concept of effective relativistic potentials dis- 
cussed in Paper I, which are both easy to calculate and to imple- 
ment in an existing Newtonian code, we require two criteria to be 
fulfilled for an effective relativistic potential Oak* including cor- 
rections due to rotation (referred to as Case Arot in the follow- 
ing). Firstly, any modification should only involve the changes of 
the gravitational potential, thereby guaranteeing that the struc- 
ture of the term -pV<I>A in the momentum equation remains un- 
changed and is still integrable, and also that the Newtonian hy- 
drodynamic equations remain unaltered. Secondly, in order to 
retain the proven accuracy of the effective relativistic potential 
<E>a, the new potential 4>Arot should reduce to that potential in the 
spherical limit. 

For Cases R and A the effective relativistic potentials are 
computed by replacing only the monopole $ N of the Newtonian 

gravitational field <&n with <I>tov or Otov.a, respectively. The ob- 
vious restriction of this approach is that any non-spherical (and 
in particular rotational) effects are accounted for exclusively in 
the higher multipoles of the Newtonian potential On, limiting the 
GR corrections to the monopole. In contrast, Case Arot also in- 
volves GR corrections to the multipoles by modifying the source 
term of the Poisson equation ©. 

Hence, we introduce a modified Newtonian potential <D N Alot , 

whose spherical part <t>N,Arot is (as in the other cases) substi- 
tuted by the corresponding modified spherical TOV potential 

^TOV.Arot: 



= <D 



N.Aiot 



-0> 



N.Arot 



+ o 



TOV.Arot ■ 



(13) 



The modified Newtonian potential <t>N,Arot is calculated by ap- 
plying several changes to the source term in the Poisson equa- 
tion (|2]i. As a first step, we replace the density p by pWn, where 
the angular Lorentz factor Wn = (1 - v^)~'/ 2 with v 2 n = v^+v 2 is 
restricted to the meridional and rotational velocity components 
vg and v^. With this choice, in the case of a purely angular veloc- 
ity the density equals the conserved density in the GR continuit y 
equation in the ADM 3 + 1 split (see e.g. iBanvuls et al.llT997b . 
Disregarding the radial velocity component v r in the Lorentz fac- 
tor is motivated by the Birkhoff theorem, according to which 
neither the gravitational mass nor the vacuum field of a spheri- 
cal star change if the matter is subject to radial motion only. 



A comparison of the spherical Newtonian potential <J>n with 
the (unmodified) spherical TOV potential O-rov suggests addi- 
tional GR corrections for the source term of <3>N,Arot- Neglecting 
the pressure term in Eq. (O we get 



2 <90 TO v _- (h) 



dr 



m TO v 



TOV 



and from Eq. (0 



2<90 N - 
r — — = m N . 



dr 



(14) 



(15) 



By identifying m-rov as a relativistic generalization of m^, we 
obtain an additional correction factor ^/F^ OVAlot , which approx- 
imately takes into account the influence of the relativistic en- 
thalpy h and the self-interaction of the gravitational field in GR 
through Fxov,Arot- The factor h can also be inferred by comparing 
the gravity terms V€>n and /iVv in the equations of rotating equi- 
librium configurations in Newtonian and GR gravity, respec- 
tively, where v is a component of the spacetime metric which 
corresponds to the gravitational potential <J>n in Newtonian grav- 
ity. 

Thus, for calculating <I>N,Arot we finally arrive at 



A<&N,Arot = 47Tp A rot, 

with the following expression for the source term: 
phW n 



PArot 



F 2 

1 TOV.Arot 



(16) 



(17) 



Accordingly, the effective relativistic potential <t>N,Arot can then 
be calculated in integral form as 



N.Arot 



I dr' d6' dip' r 
Jv 

and its spherical part as 
f M dr' 

<&N,Arot(r) = - —r 

Jo r' z 



2 • n , PArot 

sinfl ■ ■ 



r - r' 



Here the modified spherical Newtonian mass is defined as 



'«N,Arot(r) = An dr' r p Arot , 
Jo 

with the modified spherical density p A rot given by 
_ _ {p){h){W n ) 

PArot — _2 

rTOV.Arot 



(18) 



(19) 



(20) 



(21) 



Note that p A rot is not computed as the angular average of pArot 
itself but rather of its constituents, because it is a quantity derived 
from other hydrodynamic quantities. The angular average of the 
Lorentz factor is given by (Wn) — ((1 - v^y 112 ). 

For the modified spherical TOV potential of Case Arot we 
use an expression that is analogous to Eqs. (O and (flOt : 



— P°° dr' 

OTOv.ArotW = - I — ('WTOv.Arot + 4;rr' 3 (P)j ■ 

(h) 



(22) 



- TOV.Arot 
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This potential also contains several additions related to non- 
spherical effects. First of all, in the modified spherical TOV mass 

/WO = 4*r f dr' r' 2 (p W n > (1 + <e»f Ar ot, (23) 
Jo 

the angular averaged density (p) is substituted by (pWn) with 
the same rationale as in the source term of the Poisson equation 
( fToT i, while the definition of 



3. Numerical methods and models 

In order to comprehensively test the new effective relativistic po- 
tential in situations where centrifugal forces play an important 
role, we apply it to simulations of the collapse of a rotating stel- 
lar core to a PNS and of a rapidly rotating NS in equilibrium. In 
Sect. |3.1| we summarize the properties of the code used for these 
simulations, while in Sect. l3.2l we describe the setup of the initial 
models and the EoS for each of the two scenarios. 



I i / \2 2mxOV,Arot 
r T OV,Arot(>") = \ 1 + <V r > 



(24) 



follows those of Eqs. ([8]) and ( fT2b . The introduction of the cen- 
trifugal factor 



{pK) 

(pWn) 



(25) 



in the effective relativistic potential d22l) is suggested by com- 
paring the equations for stationary equilibria of rotating self- 
gravitating fluids in Newtonian and GR gravity. In essence, the 
centrifugal factor results in a larger gravitational acceleration 
for high rotational velocities. It mimics a relativistic effect that 
is, for instance, known from the motion of a test particle in a 
Schwarzschild backgrou nd spacetime, wh ere the effective rela- 
tivistic potential (see e.g.|Straumann 20041) 



V eS =ll-— III +-T 



(26) 



for the motion of a test particle with specific angular momentum 
I = rvn contains a factor (1 + l 2 /r 2 ), which is identical to (1 - 
v n) 1 to 01< ^ el ' Oiy 2 ). The density averaging is preferable to a 
volume averaging to account for the rotational flattening in cases 
of rapid rotation. 

As required the new effective relativistic potential <J>Aiot by 
construction reduces to <D>a in spherical symmetry. In this limit 
the multi-dimensional and spherical (modified) Newtonian po- 
tentials are identical, On = 3>n and i>N,Arot = 4>N,Arot- Thus, 
they cancel each other in Eqs. (0 and ( fT3l , respectively, while 
the modified TOV potentials ( [T0ll22l become identical as (Wn) 
vanishes. 

Note that different from Paper I, for reason of simplicity we 
omit all contributions due to neutrino effects in the effective rel- 
ativistic potentials discussed here. However, it is straightforward 
to add the appropriate terms to any of the potentials (see Paper I). 



2.5. General relativistic formulation using the conformally flat 
approximation 

To be able to directly assess the quality of the effective rela- 
tivistic potential approach in describing rapidly rotating compact 
objects, we have repeated all simulations using a GR hydrody- 
namics code. These results are referred to as Case GR. As in 
Paper I, we approximate the GR s pacetime metric by the con- 
forma l flatness condition (CFC; see Isenberg 1978; Wilson et al. 
Il996l) . whose excellent quality in the context of rotational stel- 
lar core collapse and a rotating NS has been demonstrated ex- 
tensively (see e .g. Shibata & Sekiguchi 20041 ICerda-Duran et al.l 
l200HlOttetaIi r2007a. and references therein). 



3.1. Code and grid setup 

All simulations are performed with the code CoCoNuT of 
Dimmel meier et al. l (l2002aL l2005) assuming axial symmetry and 
symmetry with respect to the equatorial plane. This code option- 
ally either uses a Newtonian (or alternatively an effective rel- 
ativistic) gravitational potential and Newtonian hydrodynamics 
(Cases N, R, A, and Arot) or solves the GR equations of fluid 
dynamics and the GR field equations in the ADM 3 + 1 split 
assuming CFC (Case GR). This procedure allows for a direct 
comparison of the different effective relativistic potentials with 
GR using one numerical code and an identical grid setup. 

The CoCoNuT code employs a metric solver b a sed o n 
spectral methods as described in Dimmelmeier et al. (2005). 
The (Newtonian or GR) hydrodynamic equations are formu- 
lated in conservation form, and are solved by high-resolution 
shock-capturing schemes with state-of-the-art Riemann solvers 
and piecewise parabolic cell reconstruction procedures on an 
Eulerian grid in spherical polar coordinates \r, 6). 

Rotating NS are modelled with an equidistant grid of 80 x 30 
zones. For the simulations of stellar core collapse the radial grid 
is logarithmically spaced to increase the effective resolution in 
the center where the PNS forms. In that case the grid consists 
of 200 x 30 zones in {r,0} with a central radial resolution of 
Ar c = 500 m for models with a simplified EoS, and of 250 x 45 
zones with Ar c = 250 m for the models computed with a micro- 
physical EoS. In both cases, the computational domain extends 
from r m in = to r max = 2000 km. Reflecting boundary condi- 
tions are imposed at the origin; otherwise no specific measures 
are taken to deal with the singularity at r — 0. A small part of the 
grid is occupied by an artificial low-density atmosphere extend- 
ing beyond the surface of the NS or stellar core. For more details 
about the grid setup and the b o undary conditions, we refe r to 
e.g. [D immelm eier et al.l d2005l) . iDimmelm eier et al] d2006l) . or 
IOttetalJ(l2007bl) . 



3.2. Equations of state, and models for rotating neutron stars 
or stellar core collapse 

As initial data for the simulations of NS we choose rotating equi- 
librium configurations which obey a polytropic EoS, 



P = Kp> 



(27) 



with K = 1.4553 x 10 s (in cgs unit£]) and y — 2. The models are 
comp uted with the self-consistent field method of Komats u et al.l 
( 1989). Their rotation law for the specific angular momentum j 
is given by 



j = A 2 (Q. C - Q), 



(28) 



where A parameterizes the degree of differential rotation (more 
differential for smaller values of A), and £2 C is the value of the 



This value is commonly used for polytropic NS models. 
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Table 1. Nomenclature of the iron core collapse models AxByGz (simulated with the hybrid EoS) and s20AxBy (simulated with 
the microphysical EoS). The length parameter A specifies the initial degree of differential rotation, (T/\W\)i is the initial rotation 
rate, and j\ is the stiffness of the subnuclear part of the hybrid EoS, respectively. We also give the angular velocity at the center 
for the initial configuration (Q c ,i) and at bounce (Q c ,b)- For models marked with an asterisk, we give the maximum angular velocity 
f^max.b instead of Q c b, because the density in the central region is extremely low due to the toroidal density stratification, and Q c b is 
therefore no good indicator for the rotational state of the entire core. 



Model 


A [1000km] 


(771 WD. [%1 

\ / 1 ' * lyl L 'J 


fi CJ IV 1 ] 


n c>b [ioV] 


Ti 


A1B1G1 


50.0 


0.25 


2.3 


2.9 


1.325 


A1B2G1 


50.0 


0.5 


3.1 


3.6 


1.325 


A1B3G1 


50.0 


0.9 


4.1 


6.2 


1.325 


A1B3G2 


50.0 


0.9 


4.1 


4.8 


1.320 


A1B3G3 


50.0 


0.9 


4.1 


2.8 


1.310 


A2B4G1 


1.0 


1.8 


8.0 


1.8 


1.325 


A3B1G1 


0.5 


0.25 


4.6 


3.9 


1.325 


A3B2G4 


0.5 


0.5 


6.3 


6.6 


1.300 


A3B4G2 


0.5 


1.8 


12.4 


2.3 


1.320 


A4B1G1 


0.1 


0.25 


31.1 


21 


1.325 


A4B1G2 


0.1 


0.25 


31.1 


21 


1.320 


A4B2G2 


0.1 


0.5 


43.6 


22 


1.320 


A4B2G3 


0.1 


0.5 


43.6 


7.5 


1.310 


A4B5G4 


0.1 


4.0 


139.4 


2.4* 


1.300 


A4B5G5 


0.1 


4.0 


139.4 


2.5* 


1.280 


S20A1B1 


50.0 


0.25 


1.0 


1.0 




S20A2B4 


1.0 


0.5 


9.6 


3.3 





angular velocity Q, at the center. In the Newtonian limit, this ro- 
tation law reduces to 



A 2 Q C 



A 2 ■ 



(29) 



where m is the distance to the rotation axis. 

In the simulations of rotational core collapse with a simpli- 
fied EoS the initial rotating equilibrium models are computed 
with a poly tropic EoS. They approximate an iron core sup- 
ported by electron degeneracy pressure with a central density 
p c ; = If) 10 g cm 4 , and the EoS parameters K = 4.897 x 10 14 
and y = 4/3, respec tively. For the evolution of the cores a hybri d 
EoS is employed dJanka et alJll993l : lDrm melmeier et al. 2002a). 
It consists of a polytropic contribution describing the degener- 
ate electron pressure and (at supranuclear densities) the pressure 
due to repulsive nuclear forces, and a thermal contribution which 
accounts for shock heating: 



P = P p + P th , 

where 

P P = Kpf, 



: peth(Tth - 1), 



(30) 



(31) 



with yth = 1.5. To trigger the collapse the adiabatic index is 
reduced from its initial value 4/3 to j\ < 4/3. At nuclear density 
Pnuc = 2 x 10 14 g cirT 3 the adiabatic index is raised abruptly to 
a value j2 = 2.5 to mimic the stiffening of the nuclear EoS, 
which results in a rebound of the core (the core bounce) and the 
formation of the PNS. The hybrid EoS provides for a smooth 
transi tion between th ese two density regimes (for more details, 
see lJankaetallll993l) . 

In our microphysically more detailed calculations of stel- 
lar core collapse we e mploy the ta bulate d non-zero temper- 
ature nuclear EoS by She n et al. I (119981) in the variant of 
Mareket al. (2005) which includes pressure contributions from 



baryons, leptons, and photons. Deleptonization by electron cap- 
tu re onto nuclei and f ree protons is implemented according 
to iLiebendorferl (120051) : During collapse the electron fraction 
Y e is parameterized as a function of density based on neu- 
trino radiation-hydrody namic simulations in spherical symme- 
try_JMarek et al. 2005) using the electron capture rates of 
(Langanke & Martmez-Pinedo 2000). After core bounce, Y e is 
only passively ad vected, and any furth er lepton loss is neglected. 
Again following Liebendorfer (2005), above the trapping den- 
sity ptrap = 2.0 x 10 12 g cirT 3 contributions due to neutrino radi- 
ation pressure are taken into account. As initial models we take 
the non-rotating 20 M H solar-metallicity progenitor s20.0 from 
IWooslev et al.l (120021) and impose the rotation law j28l . Note 
that these initial models are not in equilibrium, because the pro- 
genitor model s20.0 already shows considerable negative radial 
velocities and the models are not relaxed with respect to the ro- 
tation added. 

In Table [TJ we summarize the parameters of all investigated 
stellar core collapse models. They are characterized by their ini- 
tial degree of differential rotation A and by their initial rotation 
parameter (T/\W\)i, which is the ratio of rotational energy to (the 
absolute value of the) gravitational binding energy. The corre- 
sponding angular velocity Q c> j of the progenitor model is also 
shown, as well as the angular velocity £2 c b at bounce. For models 
evolved using the hybrid EoS, the stiffness j\ of the subnuclear 
EoS is also specified. 

Our selection of models includes fast-rotating cores which 
end up as PNSs with spin periods in the milisecond or sub- 
milisecond range. Lower initial rotation rates - as predicted by 
recent stellar evolu tion models incorporating magnetic braking 
(Heg er et al.ll2005l) - are required to match the spin periods of 
observed young pulars, which are typically larger than 10 ms. 
However, since our goal is to test the effective potential ap- 
proach for such rotation rates as required by the scenarios (col- 
lapsars, etc.) mentioned in the introduction, we also consider 
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Table 2. Maximum density p max ,b at bounce and final maximum density p max ,f after ring-down in units of 10 14 g cirT 3 for all 
investigated core collapse models in Case GR, or evolved with the effective relativistic potential <1>r, and <&ai i, respectively. For 
each potential, the relative deviation from Case GR (in percent) is given in parentheses, and the value closest to Case GR is indicated 
by an asterisk. Note that some centrifugally bouncing models do not reach stable values for p max ,f until very long after core bounce. 



















Pmax,f 




Model 


Case GR 


CaseR 


Case A 


Case Arot 


Case GR 


Case R 


Case A 


Case Arot 


A1B1G1 


4.71 


4.55 (-3)* 


4.32 (- 


-8) 


4.42 (-6) 


3.30 


3.53 (+7) 


3.36 (+2)* 


3.40 (+3) 


A1B2G1 


4.37 


4.26 (-3)* 


3.98 (- 


-Q) 
J > 


4.16 (-5) 


3.16 


3.32 (+5) 


3.19 (+1)* 


3.25 (+3) 


A1B3G1 


3.86 


3.78 (-2) 


3.56 (- 


-1) 


3.81 (-1)* 


2.84 


3.00 (+6) 


2.82 (-1)* 


2.97 (+5) 


A1B3G2 


4.02 


3.88 (-3) 


3.70 (- 


-8) 


3.90 (-3)* 


2.78 


2.88 (+4) 


2.78 (±0)* 


2.84 (+2) 


A1B3G3 


4.19 


4.04 (-4)* 


3.93 (- 


-6) 


4.01 (-4) 


2.99 


2.11 (+4) 


3.03 (+1)* 


3.10 (+4) 


A2B4G1 


0.64 


0.44 (-31) 


0.30 (- 


-53) 


0.58 (-9)* 


— 


— 


— 


— 


A3B1G1 


4.46 


4.27 (-4)* 


4.05 (- 


-9) 


4.25 (-5) 


3.16 


3.33 (+5) 


3.18 (+1)* 


3.26 (+4) 


A3B2G4 


3.98 


3.83 (-4) 


3.76 (- 


-6) 


3.87 (-3)* 


2.54 


2.57 (+1) 


2.54 (+0)* 


2.60 (+2) 


A3B4G2 


0.59 


0.41 (-31) 


0.33 (- 


-44) 


0.58 (-1)* 










A4B1G1 


4.64 


4.16 (-7) 


3.94 (- 


-15) 


4.30 (-10)* 


2.98 


2.96 (-1)* 


2.85 (+4) 


3.09 (+4) 


A4B1G2 


4.46 


4.12 (-8) 


3.88 (- 


-7) 


4.23 (-5)* 


2.84 


2.83 (±0)* 


2.75 (-3) 


2.96 (+4) 


A4B2G2 


4.31 


3.59 (-17) 


3.47 (- 


-20) 


4.08 (-2)* 


2.56 


2.37 (-7) 


2.30 (-10) 


2.72 (+6)* 


A4B2G3 


4.00 


4.46 (+11) 


3.38 (- 


-16) 


3.91 (-5)* 


2.30 


2.14 (-7) 


1.55 (-33) 


2.42 (+5)* 


A4B5G4 


0.47 


0.36 (-25) 


0.35 (- 


-27) 


0.51 (+7)* 










A4B5G5 


1.72 


1.26 (-27) 


1.26 (- 


-27) 


1.72 (±0)* 










S20A1B1 


3.21 


3.14 (-2)* 


3.06 (- 


-5) 


3.07 (-4) 


2.69 


2.78 (+3) 


2.73 (+1)* 


2.74 (+2) 


S20A2B4 


1.61 


1.18 (-27) 


0.89 (- 


-45) 


1.64 (+2)* 


0.94 


0.75 (-20)* 


0.62 (-34) 


1.18 (+26) 



rapidly rotating models. While rare in nature, there are viable 
formation channels for such progenitors: Even when magnetic 
torques are included, it is still possible to obtain fast rotating 
cores for massive Wolf-Rayet stars for sufficiently low mass loss 
rates dWooslev & H eger 2006) (which may be realistic for low- 
metallicity stars), although few supernovae are expected to orig- 
inate from such progenitors because of the declin ing initial mass 
funct i on. Binary evolution e ffects, i.e. ac cretion (ICantiello et alj 
120071: lYoon & Langerll2005l) or mergers dFrver & Hegerll2005l) . 
may also lead to rapidly spinning cores. Again, only a small frac- 
tion of core collapse events can be explained this way: For the 
particular case of accretion-induced co llapse, the most optimistic 
value quoted bv lDessart et alj d2007al) is only 1 .5 x 10~ 4 galactic 
events per year. 

4. Results 

4.1. Dynamics of rotational stellar core collapse 

Previous simulations, considering a large variety of rotation rates 
and profiles in the progenitor core but simplifying the com- 
plex microphysics and/or the influence of GR, found qualita- 
tively and quantitati vely different types of GW burst signals (see 
e.g. the work by Monchmever e taT] I199U IZwerger & Miiller 
1997: iDimmelmeier et al.ll200 2b). These can be classified de- 
pending on the collapse dynamics: Type I signals are emitted 
when the collapse of the homologously contracting inner core is 
not strongly influenced by rotation, but stopped by a pressure- 
dominated bounce due to the stiffening of the EoS above nuclear 
matter density. This leads to the formation of the PNS with a 
maximum core density p max > p nuc after a few ring-down oscilla- 
tions. Type II signals occur when centrifugal forces, which grow 
during contraction due to angular momentum conservation, are 
sufficiently strong to halt the collapse, resulting in consecutive 
(typically multiple) centrifugal bounces with intermediate co- 
herent re-expansion of the inner core, seen as density drops by 



often more than an order of magnitude; thus here p max < p nuc 
after bounce. Type III signals result from a pressure-dominated 
bounce when the inner core has a very small mass at bounce due 
to a soft subnuclear EoS or very efficient electron capture. 

Accordingly, we discuss the quality of the new effective rel- 
ativistic potential by first applying it to rotating stellar cores 
which collapse to form a PNS, focusing on three different cases: 
a regular pressure-dominated bounce model with slow or mod- 
erate rotation, a centrifugal single bounce model with rapid and 
strongly differential rotation (characterized by a marked toroidal 
density stratificatiorQ that develops in the core), and a multi- 
ple centrifugal bounce model. As the Type III bounce models 
have very similar dynamics to Type I models, we refrain from 
discussing them separately. The reader should refer to Table [2] 
throughout the discussion, as it presents key quantities (namely 
the maximum density at bounce and the final maximum density) 
for our entire set of models, thus providing further illustrative 
examples for moderately rotating cores with weak (Al mod- 
els, S20A1B1) and strong (most A3 and A4 models) differen- 
tial rotation, as well as multiple (A2B4G1, A3B4G2) and single 
(A4B5G5, A4B5G5 and s20A2B4) centrifugal bounce models. 

4.1.1. Slow and moderate rotation with pressure-dominated 
bounce 

For the regular collapse models with slow or moderate rotation 
we essentially confirm the findings of Paper I, as the effective 
relativistic potentials <I>Arot and <I>a are identical in the spheri- 
cal limit. For the models A1B3G3 and A4B1G2 presented in 
Paper I, as well as for all our other models with the simple hybrid 
EoS which undergo pressure-dominated bounce, €>Arot yields re- 
sults that are in very good agreement with those of GR simu- 
lations. In particular, after bounce and ring-down p max is repro- 

2 In this paper, a configuration will be called "toroidal" if the maxi- 
mum density is reached off-center. 
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J — Relativistic (Case GR) - 

a — Newtonian with regular TOV potential (Case R) ~ 

3 Newtonian with modified TOV potential (Case A) - 

S __ Newtonian with modified TOV potential ~ 

J and rotation corrections (Case Arot) _ 



t - t b [ms] 

Fig. 1. Time evolution of the maximum density p max for the 
slowly and almost uniformly rotating model A1B1G1 for 
Case GR (solid line), Case R (dashed-dotted line), Case A (dot- 
ted line), and Case Arot (dashed line), respectively. The results 
of all Newtonian simulations with effective relativistic potentials 
differ from the GR results by only a few percent, tb is the time of 
core bounce. 
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Fig. 2. Same as Fig. Q]but for the multiple centrifugal bounce 
model A2B4G1. The effective relativistic potential <f>Amt repro- 
duces the GR result to within 10% even though this model 
is strongly influenced by centrifugal forces. The additionally 
shown purely Newtonian simulation grossly underestimates 

Pmax ■ 



duced to an accuracy of 7% or better (see Table[2]i. Note that for 
these nine models (T/\W\\ < 0.1. The maximum density p maXj b 
at bounce is usually a bit underestimated, but correct to within 
10%. Typically, Case A matches Case GR most closely in the 
post-bounce phase, followed by Case Arot (the new potential in- 
cluding rotation corrections) and Case R (the potential with pure 
TOV corrections), while Case Arot works best during bounce. 
This behaviour is illustrated for model A1B1G1 in Fig. Q] The 
agreement is even better for the microphysical model S20A1B1, 
whose initial rotation rate and profile is identical to that of model 
A1B1G1. Here the error in p max (compared to the GR case) is 
smaller than 5% at bounce (where now Case R agrees best with 
Case GR) and 3% after ring-down, as summarized in Table [2] 

4.1.2. Multiple centrifugal bounce 

The limits in applying the old effective relativistic potentials pre- 
sented in Paper I (Cases A and R) become apparent for mod- 
els which collapses slowly due to j\ being close to 4/3 and are 
strongly influenced by centrifugal forces due to a considerable 
amount of initial rotation. Such models undergo multiple cen- 
trifugal bounces where always p max < p nllc , and do not settle 
down to a quasi-equilibrium state even after several bounces, 
and the resulting waveform of the emitted gravitational radia- 
tion is of Type II. It should be emphasized that multiple bounces 
do no longer occur once deleptonization is taken into account, 
essentially because a choice of y\ very close to 4/3 does not 
constitute an adequate approximation for the dynamics of the 
collapse. However, a rather high value of j\, such as to produce 
multiple bounces, still serves as an interesting test case which 
allows us to check the corrections due to rotation included for 
Case Arot in a regime where centrifugal forces dominate the 
evolution of the core around bounce. For this reason, we con- 
sider one model (A2B4G1) which provides a good example for 
the multiple bounce scenario. 



In this model (and comparable models), p max .b at the first 
bounce is underestimated by more than 50% in Case A and by 
more than 30% in Case R (see Fig. [2] and also Fig. 13 in Paper I), 
while the use of the new potential <I>A ro t reduces this error to 
less than 10%. This is an extraordinary result, as in a collapse 
situation with multiple centrifugal bounces the evolution of p max 
sensitively depends on the balance between centrifugal forces 
and the pressure gradient for an EoS with a j\ close to 4/3. 

For rotational core collapse simulations with a simple hy- 
brid EoS, the appropriate modeling of GR effects (via a consis- 
tent relativistic hydrodynamic formulation like in Case GR or 
by the use of the new effective relativistic potential <t>Arot in an 
otherwise Newtonian framework) strongly reduces the range in 
para meter space where multiple centrifugal bounce models oc- 
cur {Dimmelmeier et al. 2002b). Moreover, independent of the 
inclusion of relativistic gravity, only models with j\ > 1.31 
will exhibit such a bounce behavior, even in a purely Newtonian 
treatment. However, as will be explained in detail in the next 
Sect. 14.1.31 the occurence of models showing multiple centrifu- 
gal bounces results from the use of a simplified EoS, and van- 
ishes if an appropriate microphysical matter model is applied. 

4.1.3. Rapid and strongly differential rotation with single 
centrifugal bounce 

If the microphysical EoS and the approximate deleptonization 
scheme described in Sect. I3.2I are employed, the collapse dy- 
namics and the GW burs t signature are ex clusively of Type I. 
Th is was first observed bvlOtt et all d 200 7a). and then confirmed 
by iDimmelmeier et all (2007) who considerably extended the 
number of investigated models and comprehensively explored 
a wide parameter space of initial rotation states. They found that 
the combination of deleptonization during core collapse and (to 
a lesser degree) the influence of GR lead to this emergence of 
a generic GW signal type. In most microphysical simulations, 
the core undergoes a regular pressure-dominated bounce, which 
results in a dynamical behavior that is well represented by the 
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Relativistic (Case GR) 

Newtonian with regular TOV potential (Case R) 

Newtonian with modified TOV potential (Case A) 

Newtonian with modified TOV potential and 
rotation corrections (Case Arot) 



Fig. 3. Same as Fig. Q~|but for the centrifugally bouncing mi- 
crophysical model s20A2B4. The effective relativistic potential 
^Arot yields very good results around core bounce, but later 
causes an overestimation of p m ax- In a purely Newtonian sim- 
ulation this model does not even start to contract. 
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Fig. 4. Same as Fig. [T]but for the rapidly and strongly differ- 
entially rotating model A4B2G3. Unlike the effective relativis- 
tic potentials Ojov and OArot, potential Oa cannot reproduce 
the correct (supranuclear) maximum density of the core after 
bounce. 



evolution of p max shown in Fig. [T] Accordingly, the new effec- 
tive relativistic potential <l>Arot as well as potentials Oa and Or 
again yield excellent results, as they also do for models based on 
a simple hybrid EoS (see model S20A1B1 in Table|2]). 

Only for models with very fast initial rotation (with an ini- 
tial valued of T/\W\ > 0.01 . . .0.02), the core bounce is not 
caused by the stiffening of the EoS but rather by centrifugal 
forces. However, in contrast to the multiple centrifugal bounce 
behavior of the (simple hybrid EoS) model A2B4G1 discussed in 
Sect. 14.1.21 with microphysics only a single bounce with subse- 
quent formation of a quasi-stationary PNS is observed, reflected 
by a Type I waveform. For such cases, where rotation plays a 
similarly crucial role for the dynamics as e.g. during the disk 
formation phase in the collapsar model, the potential <&Arot now 
clearly outperforms potentials Oa and Or (see Fig. [3). At core 
bounce we find p max ,b = 1 -64 x 10 14 g crrT 3 for Case Arot, which 
is only 2% larger than the value 1.61 x 10 14 g crrT 3 obtained in 
the GR simulation, while in Cases A and R the value of p maX; b 
is underestimated by 27%, and 45%, respectively. In the post- 
bounce phase, however, <I>Arot leads to an overestimatate of p max> b 
by roughly the same amount as Or underestimates the correct 
value. In summary, for this rapidly rotating configuration of a 
core that is stabilized against gravity by centrifugal forces rather 
than pressure gradients, the new effective relativistic potential 
with rotational corrections performs better than the old poten- 
tials at bounce and at least as good after bounce. 

We also emphasize that the microphysical model s20A2B4 
does not even start to contract in Newtonian gravity due to ro- 
tational stabilization. This case where a GR instability leads to 
the contraction of a core that is still well in the Newtonian limit, 
clearly illustrates the usefulness of the effective relativistic po- 
tential approach. Such a GR instability is also responsible for 
the collapse of supermassive stars with masses between 10 4 and 
10 8 M G , which are candidates for being progenitors of super- 



3 Of course, the precise value for the initial rotation parameter be- 
yond which a centrifugal bounce at subnuclear den sities occurs also de- 
pends on the rotation profile and the choice of j\ . Dimmelmeier et al. 
(2007) should be referred to for a detailed analysis. 



massive black holes. Our results suggest that in this scenario the 
use of an effective relativistic potential should be able to capture 
the qualitatively correct collapse dynamics, too, while a purely 
Newtonian treatment fails. 

If a simple hybrid EoS is used, the situation which comes 
closest to that of the single centrifugal bounce of the microphysi- 
cal models is when fast initial rotation is combined with a subnu- 
clear adiabatic index j\ that is significantly lower than 4/3. Then 
the collapse of the core is also stopped by only a single bounce 
(albeit with possibly very strong ring-down oscillations as seen 
for model A4B2G3 in Fig. |4] or a considerable re-expansion of 
the core by sometimes more than a factor of 10 in p max as shown 
for model A4B5G5 in Fig. |5j. The post-bounce PNS is either 
exclusively or predominantly stabilized by centrifugal forces. 
Many such models have a toroidal density stratification where 
the off-center maximum density p max exceeds the central density 
p c by more than one order of magnitude. 

In such models the failure of an effective relativistic po- 
tential without corrections due to rotation becomes particularly 
apparent when only a small (possibly off-center) region of the 
core reaches supranuclear densities and the bounce is caused by 
a combination of the stiffening of the EoS and by centrifugal 
forces, as in model A4B2G3 (see Fig. |4|. Here, the potential 
Oa is unable to produce a PNS that is sufficiently compact for 
the maximum density to remain above p nuc , resulting in a post- 
bounce configuration where p max is 33% lower than in the GR 
case. The rotational corrections incorporated into the potential 
Aro t on the other hand eliminate this weakness and lead to a 
final p max which is correct to within 6%. 

The results presented in Figs. [4] and [5] demonstrate that for 
very rapidly rotating models with a simple matter prescription, 
the new potential OA ro t is superior to Oa both at bounce (show- 
ing perfect agreement with the GR results for model A4B5G5) 
and in the post-bounce phase, while it performs at least equally 
well as Or. We again emphasize that the TOV potential Or does 
not give the correct results for slow or vanishing rotation, and is 
thus not as versatile in handling a wide variety of rotation states 
as Arot . 



2.0 



B. Miiller et al.: An effective gravitational potential for rapid rotation 

-, , , , , , 1 2.5 




Relativistic (Case GR) 

Newtonian with regular TOV potential (Case R) 

Newtonian with modified TOV potential (Case A) 

Newtonian with modified TOV potential and 
rotation corrections (Case Arot) 



10 



- 1. [ms] 



Fig. 5. Same as Fig.[T]but for the single centrifugal bounce model 
A4B5G5. Again the results obtained with the effective relativis- 
tic potential <l>Arot show excellent agreement with those found in 
GR. Note that in this model the density maximum is located off- 
center, and that the differences between the Case A and Case R 
are minimal. 



In summary, we find that unlike the other effective relativistic 
potentials, <I>Arot reproduces the GR results to within an error of 
at most 15% in p max in the entire investigated parameter space of 
rotation for both simple and microphysical models. The new ef- 
fective relativistic potential also neither produces an overly com- 
pact PNS in spherical symmetry or for slow rotation - like the 
TOV potential Or - nor does it fail for centrifugal bounces as 
does <Da, and to a smaller extent also <D R . 

4.2. Structure and deformation of proto-neutron stars 

Following Paper I, we do not confine our investigation to lo- 
cal quantities like the maximum density p max , but also study the 
spatial structure of each model by comparing radial profiles ob- 
tained in GR with simulations using the new effective relativistic 
potential with corrections due to rotation. 

All effective relativistic potentials discussed in this work as 
well as in Paper I are constructed from a solution of the TOV 
equation, which is formulated in Schwarzschild radial coordi- 
nates. Therefore, our (otherwise) Newtonian simulations presup- 
pose a specific gauge choice which is different from that used in 
the GR simulations, where the metric equations in the CFC ap- 
proximation are based on isotropic radial coordinates. As a con- 
sequence, it is necessary to apply a coordinate transformation in 
order to be able to perform a resonable comparison between both 
classes of models. 

As in Paper I we take the standard Schwarzschild line ele- 
ment in spherical symmetry, 



ds 2 = -a 2 dt 2 + 



1 



■dr 2 + r 2 dQ 2 , 



(32) 



1 - 2m I r 

and replace it by the one in isotropic radial coordinates: 

K, = "« 2 df2 + <t> 4 ( d 4 + 4 d^ 2 ) , (33) 

For details about this transformation we refer to Appendix B of 
Paper I, in particular Eq. (B.8). 



Relativistic (Case GR) 
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corrections (Case Arot) 
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Fig. 6. Radial profiles of the density p for model A4B2G3 (top 
panel; 40 ms after bounce) and model A4B5G5 (bottom panel, 
at the time of bounce) at various polar angles for Case GR (solid 
lines) and Case Arot (dashed lines), respectively. For the poten- 
tial fl>Arot we a Pply a transformation to the isotropic radial co- 
ordinate (see text for details). Both models rotate rapidly and 
strongly differentially, and model A4B5B5 bounces at subnu- 
clear densities. 



As long as the shift vector ft' is close to zero and the deviation 
of the metric functions a and <p from spherical symmetry are not 
too large, rj so then yields a good estimate for the (in general not 
spherically symmetric) GR spacetime metric. These conditions 
are very well satisfied even during the most dynamic phases of 
the collapse: The radial shift vector component typically reaches 
a maximum value of the order |/8 r | ~ 0.01 (in units of the speed of 
light c) near bounce and then decreases as the PNS settles down 
to a quasi-stationary state, while the meridional component fig 
is even smaller. The rotational component f}^ is about one or- 
der of magnitude smaller than the maximum rotational velocity 
V v , i.e. its absolute value rarely exceeds 0.01 even for the most 
rapidly rotating models. Moreover, the anisotropy of the metric 
functions a and <p is significantly smaller than the anisotropy of 
the matter fields. 

The potentials Oa and 3>Arot are identical in the non-rotating 
limit and differ only slightly for slow rotation, as the corrections 
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Fig. 7. Density contours for model A4B5G5 at bounce in full 
GR (solid lines), and for Case Arot (dashed-dotted lines). log I0 p 
is indicated by the contour labels. As in Fig.|6]we apply a trans- 
formation to isotropic radial coordinates. 



due to rotation are negligible. Thus, we concentrate our discus- 
sion on the two rapidly rotating models A4B2G3 and A4B5G5. 
Furthermore, we investigate only the innermost regions of the 
collapsed stellar core constituting the nascent PNS, which is 
in approximate equilibrium after the core bounce and the ring- 
down phase. 

Figs. [6] and [8] illustrate that the radial profiles of the den- 
sity p and the specific internal energy e in Case Arot are in very 
good agreement with those obtained in GR for both models. To 
demonstrate that this excellent matching is independent of the 
angular location, we show profiles of p at three different latitudes 
- along the polar axis, at 6 — n/4, and in the equatorial plane, 
respectively. For model A4B2G3, in both cases (Arot and GR) 
we observe after the ring-down phase a similar toroidal structure 
with a maximum density of p max « (2.3-2.4) x 10 14 g ctrT 3 at 
a radius of 5 - 6 km in the equatorial plane, and a central density 
p c close to 1.0 x 10 14 g cm -3 . Even for model A4B5G5, which 
develops a very extreme toroidal structure and is depicted in the 
highly dynamic phase during core bounce in the lower panel of 
Fig. [6] the new potential <l>A ro t yields results which closely match 
those in GR. The maximum density p max is reached at r ~ 20 km 
in the equatorial plane. Fig. [7] shows density contours for this 
model down to densities of 10 10 g cm 4 , and demonstrates how 
closely Case Arot reproduces the spatial structure of that model 
in full GR. The step-like features in e in the low-density regions, 
which are caused by shock heating, are also reproduced quite 
accurately, as can be seen in Fig. [8] 

Due to the very good agreement of the density and internal 
energy profiles between Cases GR and Arot already during the 
dynamical phase of the core bounce we also infer that the use of 
the transformation from Schwarzschild to isotropic radial coor- 
dinates needs not be limited to the case discussed in Paper I, i.e. 
a spherical PNS in the late post-bounce phase, but works even 
for strongly rotating and significantly aspherical configurations. 



4.3. Structure of rotating neutron stars 

The core collapse simulations presented here are limited to a 
maximum compactness of the PNS of 2Mpns/^pns ~ 0.2 in the 
early post-bounce phase, where Mpns and /?pns are the mass and 
radius of the PNS, respectively. In order to extend the assessment 
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Fig. 8. Same as Fig. [6] but for the specific internal energy e. Note 
that potential <l>Arot (dashed lines) closely reproduces the step- 
like features in e of the model in GR (solid lines). 



of the quality of a Newtonian simulation with an effective rela- 
tivistic potential to a more strongly relativistic regime, we also 
investigate equilibrium models of a rotating NS using the poly- 
tropic EoS (1271 1, These can serve as an approximation to either 
a cooling PNS which has already evolved over several 100 ms 
(and cannot be directly simulated with our code due to compu- 
tational constraints and the missing necessary microphysics) or 
a cold NS. Such models also enable us to check the reliability of 
our approach for scenarios such as accretion-induced collapse. 
Similarly, the ability to correctly handle equilibrium models of 
compact objects with the effective relativistic potential approach 
is a prerequisite for investigating instabilities or oscillations in a 
NS. 

Therefore, we also perform an analysis of various rotating 
NS models in GR and for the different effective relativistic po- 
tentials. We find that the new potential <t>Arot reproduces quite 
closely the correct GR dependence of most global quantities 
(like gravitational mass, binding energy, etc.) on e.g. the axis 
ratio of the NS, even for strong and differential rotation, and sig- 
nificant toroidal deformation. 
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Fig. 9. Dependence of the gravitational mass M g on the axis ratio pig 1Q Same as Rg bm showing ^ dependence of the grav . 



r p /r e for rotating NS models with a maximum density p, 
3.95 x 10 14 g ctrT 3 , and a differential rotation parameter A = 0.3. 
We show the results for Case GR (solid line), Case R (dashed- 
dotted line), Case A (dotted line), and Case Arot (dashed line), 
respectively. The results obtained with the effective relativistic 
potential <I>Arot including corrections due to rotation agree best 
with the ones of GR. 



itational mass M g on the maximum mass p max for rotating NS 
models with an axis ratio r p /r e = 0.4 and a differential rotation 
parameter A = 0.3. Again the results obtained with the effective 
relativistic potential i>Arot agree best with those in GR, while in 
the purely Newtonian case the values for M„ are consistently too 
high. Without rotation (thick lines) the Cases GR and R, as well 
as A and Arot are identical, respectively. 



As long as the NS is not too compact (i.e. p c < 4 x 
10 14 g crrT 3 , which are typical values for a PNS during the first 
few 100 ms after core bounce), the results for Case R and A do 
not differ strongly. This is evident from Fig.|9j where we plot the 
gravitational mass M g for a sequence of NS models with con- 
stant central density p c = 3.95 x 10 14 g cirT 3 and differential ro- 
tation parameter A = 0.3 (where A - A/ r e ), and vary the axis ra- 
tio r p /r e . Here r p and r e are the polar and equatorial radius of the 
NS, respectively. Fig.|9]shows that even for an extreme axis ratio 
f" P /r e = 0.05, the deviation of M g from its GR value is only about 
15% for the worst potential (Case A), and much better for poten- 
tial <J)Arot- The gravitational binding energy W is somewhat more 
sensitive to the inclusion of relativistic effects; we obtain a con- 
siderably larger value of \W\ - 6.01 x 10 53 erg for Case A com- 
pared to \W\ = 4.60 x 10 53 erg in full GR (+31%). However, us- 
ing the potentials 3>tov or <J>Arot we obtain \W\ = 5.31 x 10 53 erg 
(+15%) and \W\ = 5.04 x 10 53 erg (+10%), which is reasonably 
close to the GR value. We also find good agreement for quanti- 
ties like the baryonic mass, and the central angular velocity, as 
well as for the spatial structure of the NS models (e.g. polar and 
equatorial radii, density profiles). 

For more compact configurations, we find the deviations 
from the GR models to be more pronounced, as illustrated by 
Fig- OH where we show M g as a function of the maximum den- 
sity p max for A = 0.3 and an axis ratio r p /r e = 0.4. Apparently 
the error in M g reaches 0.6 M Q for potential (Da, and 0.47 M G for 
potential <t>TOv, but does not exceed 0.15 M Q for the new poten- 
tial <l>Arot with rotational corrections. 

It is also noteworthy that for fast rotation the deviations from 
GR in Cases A and R already manifest themselves well below 
the density p max « 5 x 10 14 g cirT 3 of maximum gravitational 
mass, which is the transition point to the unstable branch. For 
non-rotating models the deviations remain small even beyond 
that point, which is then located at p max m 20 x 10 14 g cirT 3 
(see Fig. [TUJ and cp. Paper I). On the other hand, the purely 



Newtonian results overestimate M g by a factor of 2 already for 
Pmax = 8 x 10 14 g cm -3 , regardless of the strength of rotation. 
Thus, the effective relativistic potentials, and particularly the 
new one Ar ot give a clear improvement over a purely Newtonian 
treatment. As for the core collapse simulations, Case Arot is very 
robust and accurate over the entire range of rotation states with 
varying A and r p / r e . It exhibits a maximum error in M g of * 10% 
compared to more than 20% and 30% for Case R and Case A, 
respectively. 

4.4. Oscillations of neutron stars 

The effective relativistic potential approach obviously approxi- 
mates very well the collapse dynamics of rapidly rotating stellar 
cores as compared to a simulation in GR. Furthermore, it also 
gives accurate results for the structure of a nascent PNS after 
core collapse as well as of a rotating cold NS. We now subject 
this approach to an even more stringent test for its applicability 
to the post-bounce evolution of a PNS or to instability problems 
in a rotating NS. To this end, we consider the dynamics of a NS 
near equilibrium and oscillations in the linear regime. For con- 
venience, we restrict the study to non-rotating models (where 
Case A and Arot are identical), as the results discussed here are 
qualitatively independent of rotation. 

We first evolve a very compact polytropic y — 2 NS model 
with p c = 7.9 x 10 14 g cirT 3 , which is a density value typically 
reached around 1 s after core bounce. When applying a small ra- 
dial velocity perturbation to excite radial normal modes, we find 
significantly different frequencies of the eigenmodes as com- 
pared to GR if we use any of the effective relativistic potentials 
discussed here. For instance, all these potentials overestimate the 
frequency of the fundamental radial F-mode by about 45%. This 
can be clearly seen in Fig.QTl where we present the power spec- 
trum of the time evolution of p max , in which the eigenmodes can 
be identified as distinct peaks. 
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For a less compact configuration with p c = 3.95 x 
10 14 g cirT 3 , which is obtained a few 100 ms after bounce, the 
error is somewhat smaller (~ 22%) but still much larger than 
what could be expected given the excellent agreement in the 
structure of the NS model between the GR case and the use of an 
effective relativistic potential. In addition, the higher harmonics 
of the F-mode like the //i-mode or the //2-mode also exhibit this 
apparent discrepancy (see Fig.fTTTi. 

The failure of the effective relativistic potential method in 
the context of NS oscillations can be understood by comparing 
the equations of hydrodynamics in GR and the Newtonian for- 
mulatioifl Limiting ourselves to the most straightforward case 
of the TOV potential in Case R and to spherical symmetry, we 
find the following evolution equation for the radial velocity: 



~dt 



1 dP 
p dr 



3<D 



TOV 



dr 



(34) 



An analogous equation can be obtained in GR, starting from the 
formulation o f Ivan Riperld 19791) . which uses Schwarzschild co- 
ordinates like in the case of the TOV potential. After rearranging 
some terms to allow direct identifications with Case R, one ar- 
rives at 



d(a~ l v r ) 



dt 



\_dP_ _ 50tov'| 
p dr dr J 



(35) 



Unlike the Newtonian equation d34l l. the GR equation d35l l con- 
tains an additional factor given by the ratio of the lapse function 
a and the specific relativistic enthalpy h. For a stationary state, 
the left hand sides vanish and both equations become equiva- 
lent. This explains why potential cDtov ( or other modified poten- 
tials like cDa and C>Arot) can ^ e successfully applied to model the 
structure of PNS or NS models. 

On the other hand, in the general dynamic case, the ad- 
ditional factor cannot be discarded, and will, for example, be 
present in the linearized perturbation equations. Since the lapse 
function a is significantly below 1 in compact NS (with a central 
value of a c m 0.7), this factor has a considerable impact on the 
linearized perturbation equations and consequently on the eigen- 
frequency spectrum. For special cases, such as spherical pulsa- 
tions, it may still be possible to obtain simple scaling relations 
for some of the effective relativistic potentials, but this is beyond 
the scope of this work. At present, however, we conclude that the 
effective relativistic potential approach is much less suited for 
studying the frequencies of NS oscillations than for approximat- 
ing e.g. the structure or global quantities like mass and radius. 
It is conceivable that a pseudo-Newtonian treatment should be 
used with caution whenever NS or PNS oscillations play a cru- 
cial dynamica l role, e.g. in the ac oustic explosion mechanism 
suggested by Burrows et al. (2006), whose potential is still dis- 
puted dWeinberg & Ouat aert 2008). A fortiori, this also applies 
for a purely Newtonian approach. A fully relativistic treatment 
of the hydrodynamics would probably be desirable in such a con- 
text to obtain firm quantitative predictions, even though it may 
not be a necessary ingredient for a qualitatively correct model. 

4.5. Improvement of the gravitational wave extraction 

After having investigated and compared the quality of the new 
potential <3>Amt, we now focus on the issue of GW emission 



4 We remind the reader here that except for changes in the gravita- 
tional potential source terms, in the effective relativistic potential ap- 
proach the equations of Newtonian hydrodynamics remain unaltered. 



Relativistic (Case GR) 

Newtonian with modified 
TOV potential and rotation 
corrections (Case Arot) 




Fig. 11. Power spectrum of the time evolution of the maximum 
density p max for a very compact non-rotating NS model for 
Case GR (solid line) and Case Arot (dashed line), respectively. 
The f-mode and its harmonics Hi and H2 are clearly visible. 
For all investigated modes the frequencies obtained with the po- 
tential A iot are significantly higher than those in GR. 



from simulations of rotational stellar core collapse using the 
effective relativistic potential approach. In our simulations, we 
compute the radiation quadrupole moment A^g accor ding to the 
first time-integrated Newtonian quadrupole formula dFinnf l989: 
iBlanchet et ail 19901) . 



A E2 _ 



32n 3/2 d 
VT5 df 



r 3 sin6»d6»dr- 



(36) 



p [v r (3 cos 8 - 1) - 3vg sin 8 cos 8]. 



We restrict the integration to those grid cells separated from the 
initial low-density atmosphere by at least five radial zones, thus 
suppressing numerical noise due to the reassignment of cells to 
the atmosphere. The amplitude AS? is related to the dimension- 
less GW strain h in the equatorial plane at a distance r to the 
source by (see e.g. Dimmelm eier et al.l l2002b) 



n 



15 



a E2 

= 8.85x10-" ^ 10kpC 



10 3 cm r 



(37) 



Obergauli nger et all d2006l) presented a comparison of two 
core collapse models both simulated in GR and using the TOV 
potential of Case R. They found that the simulations with the 
TOV potential lead to an error of about 50% in the GW peak 
amplitude IA^Lm ( tr,at is reached around bounce) for model 
A1B3G3 despite the fact that e.g. the evolution of the p max is 
quite similar in both cases. This is confirmed by our simulations 
of a large and representative set of stellar core collapse models 
both with a simple matter model and also with more detailed 
microphysics. 

For slowly rotating models, where hydrodynamic quantities 
like p max , density profiles, etc. obtained using an effective rela- 
tivistic potential agree with those in GR very well, the errors in 
|Ajo Imax (and the entire waveform) are significant, often between 
30% and 50%. This is exemplified by model A1B1G1, whose 
waveform is shown in the top panel of Fig. Q~2] 

The quadrupole formula is derived in the Newtonian weak- 
field and low-velocity limit, and has to be generalized to a GR 
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Fig. 12. Time evolution of the GW amplitude AS? for model 
A1B1G1 extracted from a simulation in Case GR (solid), Case R 
(dashed-dotted), Case A (dotted), and Case Arot (dashed), re- 
spectively, using the time-integrated quadrupole formula. The 
top panel shows the waveforms obtained without the transforma- 
tion to the isotropic radial coordinate for the simulations with an 
effective relativistic potential. If the transformation is included, 
the difference in the GW amplitudes compared to Case GR is 
reduced considerably (lower panel). 



spacetime before a better matching to GR results can be ex- 
pected. This requires a (non-unique) choice of the radial coordi- 
nate. In the CFC metric of Case GR, isotropic radial coordinates 
are assumed. On the other hand, as already detailed in Sect. 14.21 
the effective relativistic potentials are derived from the TOV so- 
lution in Schwarzschild radial coordinates. The analogy to the 
discussion of the spatial structure of the PNS suggests that con- 
sistency between GR and Newtonian-based simulations requires 
the use of the same gauge and coordinate choice for any com- 
parison of GW amplitudes. 

We therefore apply the transformation from the 
Schwarzschild radial coordinate to the isotropic radial co- 
ordinate (see Eq. (B.8) in Paper I) when extracting GW signals 
via the quadrupole formula. This transformation effectively 
reduces the magnitude of the factor r 3 in Eq. d36l l and thus 
leads to a considerably lower amplitude, while the waveform 



shape itself remains essentially unchanged. In many cases, the 
coordinate transformation significantly improves the agreement 
with the GR results: For the slow and almost uniform rotator 
A1B1G1 (see lower panel of Fig. [12b . AS? extracted from 
a simulation with any of the effective relativistic potentials 
now deviates by less than 10% from the GR simulation, as 
opposed to more than 40% without the coordinate change. For 
the microphysical model S20A1B1, which has the same initial 
density profile and rotation rate, the effect is also dramatic: The 
deviation is reduced from 29% to 1% for Case Arot. 

We observe a similar improvement for most other inves- 
tigated models with slow and weakly differential rotation, i.e. 
where T/\W\ < 0.1 during the whole evolution: Table [3] shows 
that for the potentials Or and Oak* the deviation from GR is re- 
duced in such cases (all the Al and Bl models), and there are 
only three exceptions for potential Oa- Even for fast rotators 
(see e.g. the B3, B4, and B5 models in Table [3), we observe 
no errors in the bounce signal larger than 20% (Case Arot) and 
24% (Case R), respectively, when the coordinate transformation 
is applied. For more rapid rotators, Case A is an exception as the 
transformation actually degrades the quality of the results. This 
is to be expected, as that potential often leads to a core density 
(i.e. compactness) at bounce that is significantly below the value 
obtained with potential Or (which always gives larger densities 
than potential Oa) and potential <I>Arot (which includes correc- 
tions due to rotation). 

In summary we can infer that the combination of the poten- 
tial <I>Arot and the radial coordinate transformation in the GW 
extraction based on the quadrupole formula delivers waveforms 
that are very accurate when compared to the ones from GR sim- 
ulations for a wide range of rotation rates and profiles of collaps- 
ing compact objects. 

5. Summary and conclusions 

We have compared pseudo-Newtonian simulations of rotational 
stellar core collapse with three different effective relativistic 
gravitational potentials to reference simulations with the confor- 
mal flatness approximation of general relativity. This has been 
done using either a simplified EoS or a microphysical treat- 
ment of matter. We have been able t o improve the effe ctive 
relativistic potential appro ach used by iMarek et alj d2005l) and 
Obergaulinger et al. (2006) both in the regime of slow and rapid 
rotation. Our tests also allow us to assess the applicability of the 
different old potentials as well as the new improvement in these 
two regimes. 

In the case of slow and moderate rotation, all potentials give 
good results for the dynamical evolution. This implies that the 
generic scenario of core collapse in massive stars can proba- 
bly be treated without including any rotational corrections like 
implemented in the new potential fl&Arot, b ecause the cores o f 
these stars are not expected to spin rapidly dHeger et al.l l2000). 
Moreover, the multiple bounce scenario, which is most prob- 
lematic for the rotationally uncorrected potentials Oa and Or, is 
suppressed once deleptonization is taken into account. However, 
we have demonstrated that the gravitational wave signal can 
still be significantly improved also in these cases by applying a 
simple coordinate transformation when utilizing the Newtonian 
quadrupole formula. Our new effective relativistic potential now 
yields wave amplitudes which are within 20% of the values 
found in relativistic si mulations, as opposed to er rors of more 
than 50% reported by lObergaulinger et all (120061) . In particu- 
lar, we observe excellent agreement for the most slowly rotat- 
ing models considered, i.e. the ones closest to the predictions 
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Table 3. Absolute value of the GW peak amplitude lASjImax f° r all investigated core collapse models in Case GR, or evolved with 
the effective relativistic potential of Case R, A, and Arot, respectively. For each potential, the relative deviation from Case GR (in 
percent) is given in parentheses. 



no radial coordinate transformation transformation to isotropic radial coordinate 



Model 


Case GR 


Case R 


Cs.sc A. 


Ctisc Arot 


Case R 


Case A 


Ccisc Arot 


A1B1G1 


1 128 


17SQ i'+SR'i 
i i oy \ \ ~jo ) 


UOl ( TtV/ ) 


lUi. / I 1 PTj 


1227 (+9) 


1 170 (+4) 


1 1 01 (+0) 


A1R7G1 


2153 


^057 1+4.7) 


Z. / VJ / ^TZi7 ) 


7861 f+33^ 

Z.OU1 ^TJJ j 




1 QSS (+Q\ 

1 7JJ I \S J 


IQRO f-81 
i yo\j \ o ) 


A1R3G1 


3340 


4?4? (+77\ 


3588 


HZJ7 ^T^Z. / ) 


7871 

Z. O / 1 ^ 1 ' ) 


9S95 (-??) 

ZJ7J ^ 


ZOJ O \ 1 J } 


A1B3G2 


2094 


2767 (+32) 


2508 (+20) 


2736 (+31) 


1956 (-5) 


1838 (-12) 


1938 (-15) 


A1B3G3 


778 


1104 (+42) 


1049 (+35) 


1055 (+36) 


837 (+8) 


810 (+4) 


801 (-7) 


A2B4G1 


651 


630 (-3) 


615 (-6) 


665 (+2) 


572 (-12) 


562 (-14) 


604 (-7) 


A3B1G1 


2456 


3711 (+51) 


3238 (+32) 


3372 (+37) 


2367 (-4) 


2187 (-11) 


2213 (-10) 


A3B2G4 


598 


791 (+32) 


762 (+27) 


779 (+30) 


623 (+4) 


606 (+1) 


608 (+2) 


A3B4G2 


842 


802 (-5) 


775 (-8) 


900 (+7) 


731 (-13) 


705 (-16) 


790 (-6) 


A4B1G1 


3269 


4610 (+41) 


3746 (+15) 


4384 (+34) 


2836 (-13) 


2485 (-24) 


2652 (-19) 


A4B1G2 


2745 


3593 (+31) 


3041 (+11) 


3550 (+29) 


2382 (-13) 


2103 (-23) 


2284 (-17) 


A4B2G2 


4490 


5067 (+13) 


4393 (+2) 


5704 (+27) 


3444 (-23) 


3084 (-31) 


3581 (-20) 


A4B2G3 


3093 


3546 (+15) 


3032 (-2) 


3823 (+24) 


2597 (-16) 


2407 (-22) 


2606 (-16) 


A4B5G4 


2466 


2754 (+12) 


2725 (+11) 


2669 (+8) 


2505 (+2) 


2479 (+1) 


2345 (-5) 


A4B5G5 


3383 


3766 (+11) 


3744 (+11) 


3614 (+7) 


3476 (+3) 


3457 (+2) 


3258 (-4) 


S20A1B1 


69 


108 (+57) 


101 (+46) 


89 (+29) 






68 (-1) 


S20A2B4 


339 


330 (-3) 


310 (-9) 


353 (+4) 






321 (-5) 



of iHeger et ail (120001) in terms of rotation rates. The best over- 
all match with relativistic gravitational wave signals has been 
obtained with <l>Arot, which is identical to Oa in spherical sym- 
metry. As in this limit < &a is superior to t he original effective 
relativistic potential O r (Marek et al. 2005), the potential A rot 
is already a very attractive choice in the regime of slow rotation. 

For strong rotation we have studied more demanding test 
problems for our potentials than considered in Paper I. In or- 
der to assess in general the usefulness of the effective relativistic 
potential approach for prospective simulations of rotational in- 
stabilities in a proto-neutron star or a cold neutron star, or for the 
accretion-induced collapse scenario, we have investigated mod- 
els of a rotating neutron star, and rotational stellar core collapse 
with high rotation rates and an extremely differential angular 
momentum distribution. We have included very rapidly rotating 
core collapse models for which T/\W\ ~ 0.4 is reached, and have 
also discussed equilibrium models of a neutron star with strong 
differential rotation and axis ratios down to r p /r e = 0.05. The ro- 
tational parameters of these models were chosen such as to also 
provide a reasonable description of the accretion-induced col- 
lapse scenario, which has however not been directly simulated 
in this work. 

As in Paper I we confirm that the accuracy of Cases R and A 
declines for rapidly rotating configurations. For instance, these 
potentials lead to a gross underestimation of the maximum den- 
sity at bounce for core collapse models with a single centrifugal 
bounce, such as model A4B5G5, where the deviation to the gen- 
eral relativistic results is as large as 27%. These 'old' potentials 
also cannot reproduce the structure of a rapidly rotating neutron 
star, although they still provide a clear improvement over the 
purely Newtonian treatment. By including rotational corrections 
in our new potential <l>A ro t, we have been able to obtain much 
closer agreement with the general relativistic results in many 
core collapse simulations as well as for rapidly rotating neutron 
star models in equilibrium. In Case Arot the maximum density 
at bounce, for example, now always lies within 10% of the cor- 



rect value, and radial profiles of strongly rotating and toroidally 
deformed models also agree very well with the ones from gen- 
eral relativistic simulations. Thus, in contrast to potentials Or 
and (Da, we infer that for all practical purposes the new potential 
^Aiot works equally well from the non-rotating limit up to very 
rapidly and differentially rotating configurations. 

However, a computation of neutron star eigenfrequencies has 
revealed large errors of up to 50% for the frequencies of the F- 
mode and its harmonics, which implies that even in spherical 
symmetry some of the dynamic properties are not captured cor- 
rectly by any of the effective relativistic potentials investigated. 
These shortcomings in accurately approximating the correct pul- 
sation frequencies reveals an inherent flaw of the effective rela- 
tivistic potential approach. 

Our findings have important implications for the range of 
applicability of an effective relativistic potential in an otherwise 
Newtonian code. The rotation rate of our most rapidly rotating 
models of (proto-)neutron stars exceed values of T/\W\ » 0.3, 
which is the expect ed extreme upper limit for rotational super- 
nova core collapse (lOtt et al.l 2007b; Dimmelm eier et all [2008) 
or accretion-induced collapse (iDessart et al.ll2006l) . and lies well 
beyond the Kepler limit for uniformly rotating neutron starfl- 
While the potentials Or and Oa, which lack rotational correc- 
tions, are already quite unreliable in this regime, such high ro- 
tation rates pose no problem for the new potential Oak*. Thus, 
our new approach is clearly the preferred choice for modeling 
rapidly rotating compact objects in the framework of pseudo- 
Newtonian simulations. 
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